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We reconsider constraints on CP-violating phases in mSUGRA, sometimes also referred 
to as the Constrained Minimal Supersymmetric Standard Model. We include the recent 
calculations of Ibrahim and Nath on the chromoelectric and purely gluonic contributions to 
the quark electric dipole moment and combine cosmological limits on gaugino masses with 
experimental bounds on the neutron (and electron) electric dipole moments. The constraint 
on the phase of the Higgs mixing mass /x, |^|, is dependent on the value of the trilinear mass 
parameter, A, in the model and on tan/3. For values of \A\ < 300 GeV at the GUT scale, we 
find \0^\/tt < 0.05, while for \A\ < 1500 GeV, \9^\/n < 0.3. Thus, we find that in principle, 
large CP violating phases are compatible with the bounds on the electric dipole moments 
of the neutron and electron, as well as remaining compatible with the cosmological upper 
bound on the relic density of neutralinos. The other CP-violating phase 9a is essentially 
unconstrained. 
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Abstract 




There has been considerable progress recently in establishing strong constraints on the 
super symmetric parameter space from recent runs at LEP Jj], 0. These constraints provide 
a lower bound to the neutralino mass of ~ 40 GeV, when in addition to the bounds from 
experimental searches for charginos, associated neutralino production and Higgs bosons, 
constraints coming from cosmology and theoretical simplifications concerning the input scalar 
masses in the theory are invoked. Indeed, when all soft supersymmetry breaking scalar 
masses are set equal to each other at the GUT scale (including the soft Higgs masses), an 
ansatz typically referred to as mSUGRA, this full set of constraints excludes the low values 
of tan/3 < 2 for /i < and < 1.65 for // > 0. MSUGRA is interesting because it is simple, 
predictive and naturally provides a stable dark matter candidate |§ (an LSP bino-type 
neutralino) over most of its parameter space. 

MSUGRA (as well as the MSSM, which does not necessarily make any a priori assumption 
concerning the scalar masses other than the constraint that they should not lead to charge 
and color breaking minima @, 0) is well known to contain several independent CP- violating 
phases. If we assume that all of the supersymmetry breaking trilinear mass terms, are 
equal at the scale Mx at which the mass parameters unify, then the number of independent 
phases reduces to 2, which we can take as 9 a and 9^. In principle, the supersymmetry 
breaking bilinear mass term B can also have a complex phase. However, the phase of \i 
can always be adjusted so that 9b = — 9 ^ by rotating the Higgs fields so that their vacuum 
expectation values are real |J, and so there are only two independent physical phases, which 
we take as above. It is well known that these phases can lead to sizable contributions to the 
neutron and electron dipole moments 0. To suppress the electric dipole moments, either 
large scalar masses (approaching 1 TeV) or small angles (of order 10~ 3 , when all SUSY masses 
are of order 100 GeV) are required. It has been commonplace to assume the latter, though 
the possibilities for large phases was considered in || and more recently in || . To reconcile 
large phases with small electric dipole moments, some of the sparticle masses are required 
to be heavy. In ||, either large sfermion or neutralino masses (or both) were required. 

However, unless i?-parity is broken and the LSP is not stable, one requires that the 
sfermions be heavier than the neutralinos. If they are much heavier, this would result in an 
excessive relic density of neutralinos. In fTD"! , we showed that there was a strong correlation 
between the CP-violating phases in the MSSM, the cosmological LSP relic density and the 
neutron and electron electric dipole moments, when all of the parameters are set at the weak 
scale and hence are independent of any RGE evolution. Large CP violating phases can also 
have an important impact on elastic scattering cross-sections used to calculate dark matter 
detection rates ||11|| . In |P2 |, we further showed that in mSUGRA, by combining cosmological 
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constraints on the mass of the LSP with experimental bounds on the neutron and electron 
EDM's, we can bound one of the two new phases in mSUGRA to lie within < n/10 for 
Aq = \A(Mx)\ ^ 1 TeV. The calculation of the contributions to the neutron and electron 
electric dipole moments has recently been reinvestigated in || [TB|] correcting some errors 
in H and including the contributions from the chromoelectric and purely gluonic operators 
(for the neutron only) which were previously thought to be small . These operators, can 
in fact lead to cancellations which enlarge the allowable range for the CP violating phases. 
Some analytic approximations to the cancellations in the edms were done in [|HJ. Here, we 
wish to reexamine this problem in light of the recent calculations of [§, 13 1 and to address 
some questions regarding the independence of the two phases when the RGE evolution is 



taken into account 16 



In mSUGRA, once the gaugino, soft scalar masses, A and 5-terms and phases are given 
at Mxi they can be RGE evolved to the electroweak scale. In practice, we use the one- 



loop RGEs for the masses and two- loop RGEs for the gauge and Yukawa couplings [17 



The structure of the equations for the gauge couplings, gaugino masses and the diagonal 
elements of the sfermion masses are such that they are entirely real. The evolutions of the 
Ai, however, are more complicated, as the Ai pick up both real and imaginary contributions. 
For example, the evolution of A t is given by 

^ = (- Y 9l Ms - 3gl M 2 -f g\ M x + h\A h + Qh*A t ) (1) 

As one can see, A t receives real contributions qM proportional to the gaugino mass (whose 
coefficients q are different for each sfermion in a generation) and (in principle complex) 
contributions dihjAf from the heavy generation (whose coefficients di differ for the first two 
and the third generations); the phases (and magnitudes) of the Ai must therefore be run 
separately. At one loop, the evolution equation for /i is given by 

f = ^(-3^-^ + ^ + 3/^ + 3^) (2) 

and the phase of \i does not run. As it will be important for our discussion below, we also 
give the RGE for the bilinear mass term B 

^ = ^2 (r^l M 2 - g\ Mi + h\A T + 3h 2 b A b + 3h 2 t A t ) (3) 

^From Eq. we see that the phase of A will induce a phase in B (and hence in /x, when we 
rotate the Higgs fields to ensure their vacuum expectation values are real). For this reason, 
it was argued [16| that a tight constraint on 9^, required either a similarly tight bound on 
9a or a fine-tuning of the GUT value of 9b- 
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To be more explicit, the Higgs superfield if 2 (which gives mass to up- type fermions) can 
be phase rotated in such a way so as to insure real expectation values for the Higgs scalars. 
The rotation changes the phase of if 2 by an amount — {9^ + 9 b)- Not only is the phase 
of \i now fixed at 9^ = —9b, but also the initial phase of \x is physically irrelevant as it is 
canceled by the rotation. As such, one might worry, that the running of the phase of B, 
may induce a large phase for \x when 9a is large, and can be excluded by the limits on the 
electron and neutron edms [|l|J. However, since 9^ = —9b, the value of 9^ at the weak scale 
will depend on both 9b and 9a at the GUT scale. Thus, any value of 9^ may be obtained 
with a suitable choice of B at the GUT scale for any 9a- Indeed, for large 9a, 9b at the GUT 
scale is sometimes required to be comparable to 9^, is sometimes considerably larger, but is 
rarely required to be set to a precision of better than 20% to obtain an acceptable value of 
9 p. We will return to this point below with some explicit results. 



Previously [10], it has been shown that the presence of new CP- violating phases may 
have a significant effect on the relic density of bino-type neutralinos. The dominant channel 
for bino annihilation is into fermion anti-fermion pairs. However, this process exhibits p- 
wave suppression, so that the zero temperature piece of the thermally averaged annihilation 
cross-section (which is relevant for the annihilation of cold binos) is suppressed by a factor of 
the final state fermion mass 2 . This significantly reduces the annihilation rate and increases 
the neutralino relic density. Mixing between left and right sfermions lifts this suppression 
to some extent by allowing an s-wave contribution to the annihilation cross-section which 
is proportional to the bino mass 2 , but the presence of complex phases in the off-diagonal 
components of the sfermion mass matrices dramatically enhances this effect. 



(4) 



To fix our notation, we take the general form of the sfermion mass 2 matrix to be [IS 

,Ml + m)+COS 23 ( T 3l -Q^e^Ml - m/ ^ \ 

V -m f m s e'^f M 2 R + m) + cos 2(3Q f sin 2 9 W M\ J 

where Mum are the soft supersymmetry breaking sfermion mass which we have assumed 
are generation independent and generation diagonal and hence real. Due to our choice of 
phases, there is a non-trivial phase associated with the off-diagonal entries, which we denote 
by m f (m f e nf ), of the sfermion mass 2 matrix, and 

m f e*" = RffJt + A} = Rf\n\e ie » + \A f \e~ i(>A f , (5) 

where mj is the mass of the fermion / and Rf = cot/3 (tan/?) for weak isospin +1/2 (- 
1/2) fermions. We also define the sfermion mixing angle 9f by the unitary matrix U which 
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diagonalizes the sfermion mass 2 matrix, 

/ cos Of sin 9t e %ls \ 

U=( . )■ (6) 

V — sin Of e 11 'f cos Of ) 

The effect of CP-violating phases on the annihilation cross section is manifest through 
both the sfermion mixing angle and a combination of the CP-violating phases 

sin 2 (20/) sin 2 7/ (7) 

For the sfermion partners of the light fermions, 

2 2 

8111 i - 26f) ~ (M 2 L - M 2 R + 2Q /C os2/W^M 2 ) 2 (8) 

Therefore, while in the more general MSSM where there is more freedom to choose Ml and 
Mr, the effect of the phases on the relic density can be quite significant, in mSUGRA, the 
effect is substantially reduced since sin 2 (2$/) is small for the lighter of generations. (For 
an LSP with a mass less than m t , the effect of the phases in the stop mass matrices is 
unimportant.) 

In addition, because of fixed relationship between the many mass parameters in mSUGRA, 
the cosmological relic density is able to furnish a strong constraint on the remaining SUSY 
parameter space. For example, in Figure 1, we show the m^-min parameter plane for 
tan/3 = 2, the regions for which the relic density takes the values 0.1 < Qh 2 < 0.3. The 
upper bound comes from the requirements that the age of the universe tjj > 12 Gyr and that 
Q < 1. The choppy region at lower values of rrixn shows the effects of the Higgs and Z poles 
on the annihilation cross-section. Also shown on the figure are the current LEP2 slepton 



mass bound]19|, which is about 84GeV for large mf — m^, and a chargino mass contour 
of 91 GeV, which approximates the LEP2 bound on the m 1( / 2 — m plane except near the 
intersection of these two contours (see Jl| for detail). The hatched regions in Figure 1 is ruled 
out because it leads to a stau as the LSP. As one can see, the requirement that flh 2 < 0.3 
gives us a constraint on both m and mi/ 2 : m < 150 GeV, and mu2 ^ 450 GeV. The shape 
of the allowed region is insensitive to tan/3 for small to moderate tan/?, and these bounds 
on mo and m\/2 do not vary up to tan/3 ~ 8. At tan/3 = 10, the bound on mo is relaxed to 
~ 170 GeV, and between tan/3 ~ 15 and 20, s-channel pseudo-scalar annihilation becomes 
important, and the bounds on both m and are significantly relaxed. However, the 
EDMs are typically very large at these large values of tan (3 and yield correspondingly tiny 
upper bounds on 6* M . The cosmological region is also insensitive to A , which we take equal 
to in this plot; the region with m f < will vary somewhat with A at large tan/3, but 
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this has only a small effect on the upper bound on my 2- In Fig. |1] we've taken 9^ = 0, but 
very similar bounds on mo and my^ apply for other values of 9^. In the shaded regions, the 
lightest neutralino is mostly bino||, and for large values of the unified gaugino mass my 2, 
Xi is almost pure bino. The cosmological constraints from this figure will serve as a basis 
for our constraints on the CP-violating phases. 

The electric dipole moments of the neutron and the electron has contributions coming 
from neutralino, chargino, and gluino exchange to the quark electric dipole moment. These 
were calculated in detail in || and more recently in 13 1 correcting some sign errors. As 
we will see, this sign change plays an important role in determining cancellations among 
the different contributions to the edms. The necessary CP violation in these contributions 
comes from either in the sfermion mass matrices or 9^ in the neutralino and chargino mass 
matrices. Full expressions for the chargino, neutralino and gluino exchange contributions are 
found in [§, ||. The neutron edm has additional contributions to its edm coming from the 
quark chromoelectric operator, and a purely gluonic operator. For all SUSY mass scales of 
0(100) GeV, it was shown in [2(| that these latter two contributions to the quark electric 
dipole moment are small, coming in with the ratio 7 : O qc : Oq = 21 : 4.5 : 1 and hence 
they were neglected in our previous analysis. However as has been argued by Ibrahim and 
Nath ||, [TI|, these ratios are not general over the interesting supersymmetric parameter 
space, and in fact, the latter two contributions can in some instances even dominate. 

The contributions to the quark electric dipole moments from the individual gaugino ex- 
change diagrams are proportional to m^/m- where mi is the mass of the neutralino, chargino 
or gluino exchanged. The edm falls as m\/2 is increased, because the squark masses 2 re- 
ceive large contributions proportional to m\j 2 during their RGE evolution from Mx to Mz- 
Roughly, 

m \ ~ m l + 6m i/ 2 + C(M|). (9) 
The chromoelectric contribution has a similar dependence on SUSY masses. The gluonic 
contribution goes as 1/mJ but still falls as m%/2 is increased if gaugino unification is assumed. 
Thus even for large values of the CP violating phases, one can always turn off the quark 
electric dipole moment contributions to the neutron EDM (and similarly turn off the electron 
EDM) by making my 2 sufficiently large 0; however one must still satisfy the cosmological 
bounds discussed above. Experimental bounds are \d n \ < 1.1 x 10 _25 ecm for the neutron 
electric dipole moment and \d e \ < 4 x 10~ 27 ecm [22| for the electron EDM. Note also that the 
squark masses m~ are only weakly dependent on mo in the cosmologically allowed regions of 
Figures 1 and 2, particularly in the regions of large my 2, and so the quark EDM's will also 
be independent of m . 
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We find that, taken independently, the current experimental bounds on the electron EDM 
are generally more restrictive than the current bounds on the neutron EDM. We compute the 
electron EDM in mSUGRA as a function of 9^, 9a, and m\/2 for fixed Aq, mo, and tan f3. We 
find that both the neutralino and chargino exchange diagrams can contribute significantly to 
the electron EDM, and in fact significant cancellations between the two contributions allow 
for much larger values of 9^ than would be permitted in the absence of cancellations. In 
Fig. |2|, we display as a function of 9^ and 9 a the minimum value of required to bring 
the electron EDM below the experimental bounds for tan/5 = 2, taking mo = 100 GeV. 
We exclude points which violate the current LEP2 chargino and slepton flnj mass bounds. 
All the contour plots we present are computed on a 40x40 grid in 9^ and 6a, and features 
smaller than the grid size are not significant. Due to cancellations, the electron EDM does 
not decrease monotonically with mm in this region of parameter space, but there is still a 
minimum value of my 2 which is permitted. The cancellations also enhance the dependence 
of rriy 1 ^ on Aq. Since the chargino contribution is proportional to sin 7/, the value of 8^ 
at which these delicate cancellations occur is sensitive to A , and larger values of A cause 
the cancellations to occur at larger values of 9^. This is demonstrated in Fig. ^|(a-c), where 
contours of m^ n = 200, 300, and 450 GeV are displayed for A = 300, 1000, and 1500 GeV 
respectively. In the zone labeled "I" , < 200 GeV, while the zones labeled "II" , "III" , 

and "IV" correspond to 200 GeV < mf 2 n < 300 GeV, 300 GeV < mf 2 n < 450 GeV, and 
rriy 1 ^ > 450 GeV, respectively. For A =0, the contours are of course straight vertical, and 
the mffi = 450 GeV contours lie at 6 J it = ± 0.009. 

The extent of the bowing of the contours in Fig. ^| increases with the the fineness 
of the cancellation between the neutralino and chargino contributions, and at the larger 
values of 9^, these cancellations can be very severe, at the level of one part in 20 for O^/ir 
near 0.05 in Fig. |a, to the level of one part in a hundred or so in the most extreme cases 
in Fig. 0c. Since the two contributions scale very similarly with my2, the cancellations 
typically occur over a broad range in m^. In the extreme cases the latter is no longer true, 
as even a tiny relative shift in the magnitudes of the two contributions is sufficient to spoil 
the delicate cancellations necessary to yield a sufficiently small electron EDM at large # M 
and 9 a- Numerically, in Fig. ||a, the allowed regions are greater than 90 GeV wide in my 2 
essentially everywhere in the figure, including at large 9^ where rriy 2 n is less than 200 GeV. 
In Fig. 0b, the allowed regions in zones I and II are between 20 and 40 GeV wide in for 
essentially all 9^/tt > 0.14. And in Fig. ^c, the allowed regions in zones I and II are between 
10 and 20 GeV wide for all 9^/tt > 0.24. 

Comparing with Fig. 1 , we see that regions for which my 2 n is greater than about 450 GeV 
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are cosmologically excluded. For Fig. |2|a, for example, this restricts [Opl/ir < 0.05 for Aq = 
300 GeV (there is an identical region with 6 a — > —6a, 6^ — ► —6^), while for v4 = 1500 GeV, 
6^/n can take values as large as 0.3. These large values of 6^ are significantly different from 
the results found earlier [12]. The differences are partly due to the sign correction in the 
neutralino contribution [IB], [| to the electron EDM, which shifts the region of cancellation 
between the neutralino and chargino contributions, and partly due to the fact that here we 
consider larger values of Aq. Of course for larger values of 6^, the cancellations become finer 
and finer, and the range of over which the electric EDM bounds are satisfied become 
smaller and smaller. The neutron EDM bound will then play a role in excluding some of 
these small regions, and we will see that the regions which still remain at large Aq and 6^ 
after imposing both the neutron and electron EDM bounds have a very small extent in m^. 
Consequently, we neglect the tiny regions in my 2 which become available at larger 6^ when 
Aq > 1500 GeV. 

We can now explicitly check the degree of fine-tuning needed to have simultaneously 
a large value for 6a and an acceptable value for 6^. For the case just considered, when 
Aq = 300 GeV and 6 A = 0.5 vr, 0.02 < OJ-k < 0.05 implies that 0.03 < 6 b /tt < 0.05 at the 
GUT scale. Namely, 6b is comparable in magnitude to 6^ and must be set only to within 
40%. As we can plainly see, there is no undue fine-tuning of 6b which is necessary. For 
this value of Aq, the change in B as it is evolved from the GUT scale to the electroweak 
scale is much smaller in magnitude than B itself, so 6b cannot change much during its 
evolution, and 6b consequently must initially lie in a wedge lying in the vicinity of —6^. 
Similarly, when Aq = 1500 GeV, we see from Fig. ||c that 0.15 < 6^/tt < 0.3, which requires 
0.18 < #b/7t < 0.21. Since B runs more for larger Aq, there is actually a finer adjustment 
needed to get the appropriate 6^ in this case, even though the actual range of 6^ is larger, 
but here too the adjustment needed hardly constitutes a fine-tuning. For Aq = 1000 GeV, 
the corresponding numbers for 0.1 < 6^/n < 0.17 are 0.10 < 6b/tc < 0.14. 

We have also computed the neutron EDM in mSUGRA as a function of 6 ^,6 a, and m 1( / 2 
for fixed AQ,rriQ, and tan/3, using the naive quark model (for the effect of other choices 
for the quark contributions to the nucleon spin, see ||23|| ). In practice we find that the 
dominant contribution to the neutron electric dipole moment in mSUGRA comes typically 
from the chargino exchange contribution to the quark EDMs. However, the gluino exchange 
contribution to both the quark electric dipole and chromoelectric dipole operators can in 
some regions of parameter space yield neutron EDMs contributions of comparable size and 
opposite sign to the chargino exchange contribution, and these cancellations permit even 
larger values for albeit over a quite narrow range in m^. In Fig. |2]d, we display 
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contours of m^™ allowed by the neutron EDM constraint alone for A = 1500 GeV and 
itlq = 100 GeV. In contrast to the electron EDM, 9^ is completely unconstrained, although 
the allowed regions are very small, less than 10 GeV wide in mm for 9^/tt > 0.3 and typically 
less than 5 GeV wide for 9^/tt > 0.4. The constraints from the neutron EDM are weaker than 
the constraints from the electron EDM; however, since the cancellations in contributions to 
the electron and neutron EDMs generally occur for different ranges of m^, it is useful to 
combine the two constraints. 

Fig. |H displays contours of m?/™ f° r the same parameters as in Fig. but now requiring 
that both the electron and neutron EDM bounds be satisfied. In Fig. |3]a, the contours 
are nearly identical to Fig. ||]a, as the neutron EDM provides little additional constraint. 
Figs. |3]a and |]b, by contrast, are noticeably altered. Essentially all of the regions with 
rriy^ < 200 GeV are removed. A gap appears in the center of the allowed region in 9 a and 
6^ , and the gap grows as A is increased. Further, the width in my 2 of the allowed regions 
decreases. In Fig. |3]c, the width in myi is between 5 and 15 GeV for the entire region with 
rriy^ < 300 GeV and 9^/tt > 0.1. In Fig. [|b, the width the allowed region in zone II varies 
principally between 10 and 40 GeV. 

Thus our general conclusion when taking into account the limits from both the neutron 
and electron edm: we find the angle 9a to be unconstrained, that is, it may take any value 
between and 7r and still remain consistent with the neutron and electron edms, as well as 
with cosmology. However, values of 9a = tt/2 require somewhat specific non-zero values of 9^ 
(see Figs. 0, 0). As discussed above we do find constraints on 6> M which are quite dependent 
on the value of A Q (\A\ at the GUT scale). This can be seen by comparing Figs. ^|a- |2|c and 
[|]a- §c, which show the contours of m^j in the 9^ — 9a plane for mo = 100 GeV, tan/3 = 2. 
Recalling from Fig.|l| that < 0.3 requires mi/2 < 450 GeV, one can read off the limits 

on 9^ from the figures. In [15| it was checked that no additional constraint due to the 
CP-violating e parameter could be placed on the MSSM phases. 

The SUSY contributions to the electron and neutron EDMs also depend on tan f3. The 
contours of constant m™™ m Figs. ||, |3] keep the the same general shape as tan j3 is increased, 
but as the EDMs do tend to increase with tan/3, the scale in 9^, and the upper bound on 
6^, decreases. In Fig. |]a we display the upper bound on 9^ coming from the electron EDM, 
as a function of tan/3 for four different sets of Aq and m®. We display the same in Fig. |4]b, 
but here we require that both the neutron and electron EDM constraints be satisfied. One 
sees that varying m has only a small effect on the 9^ bounds in Fig. |]a, and a negligible 
effect in Fig. [|b. 

In summary, we have combined cosmological bounds on gaugino masses with experi- 
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mental bounds on the neutron and electron electric dipole moments to constrain the new 
CP-violating phases in the Constrained Minimal Supersymmetric Standard Model. We find 
that the inclusion of the chromoelectric and purely gluonic contributions to the quark electric 
dipole moment as calculated by Ibrahim and Nath R O, while providing significant correc- 



tions to the neutron EDM in regions of the parameter space, do not significantly affect the 
upper bound on the phases in mSUGRA, from the combination of the neutron and electron 
EDM and relic density constraints. The phase of the supersymmetric Higgs mixing mass is 
constrained by |0 m |/tt < 0.3 for A < 1500 GeV and |0 m |/tt < 0.05 for A < 300 GeV. In 
addition, there is no bound on the phase 6a of the unified trilinear scalar mass parameter A. 
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200 




Figure 1: Contours of constant VL~h 2 = 0.1 and 0.3, as a function of mo and M, for 
A = GeV and tan/3 = 2. The dotted line represents the current LEP2 slepton exclu- 
sion contour, and the dot-dashed line corresponds to a chargino mass of 91 GeV. The shaded 
region at bottom right yields a stau as the LSP. 
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Figure 2: Contours of roT, the minimum my 2 required to bring the electron EDM below 
experimental bounds, for tan/3 = 2,m = 100 GeV, and a)A = 300 GeV, h)A = 1000 GeV, 
c)A = 1500 GeV and d) as in c) but for the neutron edm. The central light zone la- 
beled "I" has mf* 1 < 200 GeV, while the zones labeled "II", "III", and "IV" correspond to 
200 GeV < m^ n < 300 GeV, 300 GeV < m^ n < 450 GeV, and m^ n > 450 GeV, respectively. 
Zone IV is therefore cosmologically excluded. 
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Figure 4: The maximum values of 6 ^ allowed by cosmology and a) the electron EDM and 
b) both the electron and neutron EDM's, as a function of tan (3, for several combinations of 
m and A . 
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